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We analyze the characteristic features of traffic breakdown. To describe this phenomenon we apply 
to the probabilistic model regarding the jam emergence as the formation of a large car cluster on 
highway. In these terms the breakdown occurs through the formation of a certain critical nucleus in 
the metastable vehicle flow, which enables us to confine ourselves to one cluster model. We assume 
that, first, the growth of the car cluster is governed by attachment of cars to the cluster whose 
rate is mainly determined by the mean headway distance between the car in the vehicle fiow and, 
may be, also by the headway distance in the cluster. Second, the cluster dissolution is determined 
by the car escape from the cluster whose rate depends on the cluster size directly. The latter is 
justified using the available experimental data for the correlation properties of the synchronized 
mode. We write the appropriate master equation converted then into the Fokker-Plank equation 
for the cluster distribution function and analyze the formation of the critical car cluster due to the 
climb over a certain potential barrier. The further cluster growth irreversibly gives rise to the jam 
formation. Numerical estimates of the obtained characteristics and the experimental data of the 
traffic breakdown are compared. In particular, we draw a conclusion that the characteristic intrinsic 
time scale of the breakdown phenomenon should be about one minute and explain the case why the 
traffic volume interval inside which traffic breakdown is observed is sufficiently wide. 

PACS numbers; 45.70.Vn,02.50.Ey,05.70.Fh,89.40.+k 
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I. INTRODUCTION. TRAFFIC BREAKDOWN 
AS THE NUCLEATION PHENOMENON 

The spontaneous formation of traffic jams on highways 
has attracted attention for the last years because of two 
reasons. The former is the importance of this problem 
for traffic engineering especially concerning the feasibil- 
ity of attaining the limit capacity of traffic networks and 
quantifying it. The latter is a fact that the jam for- 
mation nicely exemplifies the existence of various phase 
states and transitions between them in statistical systems 
comprising elements with motivated behavior, which is a 
novel branch of physics. According to the modern notion 
proposed by Kerner (see, e.g. Ref. |^, |[ ^) based on the 
experimental data |[ ^, 0] the jam formation is of suffi- 
ciently complex nature. In particular, it proceeds mainly 
through the sequence of two phase transitions: free flow 
(F) synchronized mode (S) stop-and-go pattern (J) 
[0. Both of these transitions are of the first order, i.e. 
they exhibit breakdown, hysteresis, and nucleation effects 

. The F ^ J transition can occur directly if only the 
synchronized mode is suppressed by a road heterogeneity 
Q . The recent analysis of single- vehicle data by Neubert 
et al. in particular, conformed these features and 
also discovered fundamental microscopic properties dis- 



tinguishing the synchronized mode from the other traffic 
states. 

Theoretical description of the jam formation is far from 
being developed well because the synchronized mode pos- 
sesses extremely complex structure [||. For example, it 
comprises a certain continuum of quasistable states, so, 
matches a whole two-dimensional region on the phase 
plane "vehicle density - traffic volume" in contrast to 
the free flow state. However, tackling the question of 
how to regulate traffic flow on highways, for example, by 
controlling the speed limitation in order to prevent the 
jam formation we may rough out the problem. Indeed, 
for this purpose it is sufficient to analyze the conditions 
giving rise to jams rather then the jam evolution itself. 
Such a standpoint is justified, in part, by the aforemen- 
tioned phase transitions being of the first order. The 
free flow state, presumably, should have the feasibility to 
exist at the given car density or inside its certain neigh- 
borhood. This might be the necessary requirement for 
jam formation at a fixed vehicle density, or for appear- 
ance of both the jam phase and the synchronized mode 
at a fixed traffic volume. Second, the jam formation pro- 
ceeds via the nucleation mechanism but not in a regular 
manner. Therefore, the key point in the emergence of a 
jam is the random occurrence of its critical nucleus inside 
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the synchronized mode or free flow. 

The jam formation manifests itself in the trafflc break- 
down, i.e. in a sharp drop of the traffic volume to a sub- 
stantially lower value. Detecting these events one can get 
the rate of the critical nucleus generation depending on 
the road conditions and the traffic state. In this way the 
main attention is shifted to the experimental and theo- 
retical analysis of the probabilistic features of jam forma- 
tion regarding the characteristic mean values of the trafflc 
volume as phenomenological parameters |l^, ^ 
Such a probabilistic description of the trafflc breakdown 
is the main purpose of the present paper. 

At the first glance the problem seems hopeless until 
the model of the synchronized mode is developed. Nev- 
ertheless, there are circumstances enabling us to make a 
step towards this problem right now (see also Ref. 13). 
The matter is that the F S transition is of another na- 
ture than the S ^ J transition. The former is due to a 
sharp decrease in the overtaking frequency, giving rise 
to the synchronized mode, whereas the latter is caused 
by the pinch effect (see, e.g. Ref. |l|, ||, ||). Thereby the 
main case of the F — > S transition is not fluctuations in 
the vehicle density and velocity but in other characteris- 
tic parameters of the traffic flow (cf. also Ref. By 
contrast, just these fluctuations give rise to the jam emer- 
gence in the synchronized traffic flow. As a result, the 
threshold of the F — > S transition turns out to be remark- 
ably less than that of the jam formation and attained 
at lower values of the vehicle density. So, the genera- 
tion rate of critical nuclei for the former transition has 
to be great in comparison with the latter one. Thus, on 
time scales characterizing the occurrence of the jam crit- 
ical nuclei the traffic state with respect to the transitions 
between the free flow and synchronized mode is quasista- 
tionary. Therefore the formation of a jam critical nucleus 
is the leading nonequilibrium process limiting the trafflc 
breakdown. The latter feature allows us to conflne our 
consideration solely to the jam nucleus generation and 
to regard the synchronized mode and the free flow phase 
(if they coexist in the case under consideration) as one 
trafflc state. Moreover, since a jam forms actually in- 
side the synchronized mode where the vehicle motion at 
different lane is strongly correlated we may apply to a 
single-lane road approximation that treats all the cars 
moving at different lanes on a multilane highway and be- 
ing neighboring across the highway as a single effective 
macrovehicle consisting of many cars. The macrovehi- 
cle concept is partly justifled by the empirically observed 
fact that fluctuations in the downstream flow leaving a 
freeway bottleneck can proceed without the traffic state 
change even their amplitude attains 30% of the mean 
trafflc volume |, ^ |ll|, [l2[ In any case fluctua- 

tions in the trafflc flow volume near its breakdown are 
of macroscopic nature and the critical nucleus of trafflc 
jam has to include many vehicles. This feature is also 
pointed to by the observed breakdown near a ramp-on 
occurring each time after a large vehicle cluster entered 
the freeway stream 
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FIG. 1: Illustration of the speed autocorrelation vs the num- 
ber of cars passed a fixed detector that was experimentally 
observed by Neubert, Santen, Schadschneider, and Schreck- 
enberg M. 



It should be pointed out the real structure of congested 
trafflc near a highway bottleneck is sufficiently complex, 
it contains the region of synchronized mode located in 
the close proximity of the bottleneck, the preceding up- 
stream region of moving narrow jams transformed into 
wide jams [^] . However it is quite reasonable to consider 
this structure as being induced by the trafflc breakdown 
processes arising inside the "head" of this complex jam, 
in the region of synchronized mode adjacent the bottle- 
neck. Therefore the main characteristics of the break- 
down phenomenon may be related to intrinsic processes 
taking place inside the synchronized mode on not to large 
spatial scales. The latter justifies our attempt to describe 
trafflc breakdown ignoring the complex spatial structure 
of the metastable trafflc state inside which critical jam 
nuclei originate. 

Processes similar to the traffic breakdown are widely 
met in physical systems. For example, water conden- 
sation in supersaturated vapor proceeds via formation 
of small atom clusters of critical size. Keeping in mind 
this analogy between the traffic breakdown and the phase 
transitions in physical systems Mahnke et al. jl^, [l^l pro- 
posed a kinetic approach based on the stochastic mas- 
ter equation describing the jam formation in terms of 
the attachment of individual cars to their cluster. How- 
ever, the particular form of the developed master equa- 
tion does not allow for the jam formation being the first 
order phase transition and, thus, the traffic breakdown. 
In the present paper we generalize this kinetic approach 
to describe the latter phenomenon. 

However, before passing directly to the model state- 
ment we recall the experimental data enabling us to es- 
timate the characteristic size tiq of vehicle clusters that 
are small enough so the behavior of drivers inside them 
seems to be special. From our point of view the mul- 
tilane correlations in the vehicle motion are due to the 
drivers taking into account the behavior of all the cars, 
including also the cars at the neighboring lanes, that are 
inside the region accessible to observation. Therefore the 
synchronized mode has to exhibit strong correlations in 
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FIG. 2: Schematic iUustration of the cluster transformations. 

this region. Figure |l] shows the speed autocorrelation 
function vs the number of cars passed a fixed detector 
that was experimentally found in the synchronized mode 
lH (see also Ref. |l^). We see that the car velocities are 
strongly correlated over scales spanning some ten vehi- 
cles, i.e. a car cluster of this size, tiq ^ 20, makes up 
actually a certain fundamental unit of the synchronized 
mode. In the free flow no such long-distant correlations 
have been observed. 



II. PROBABILISTIC MODEL FOR THE CAR 
AGGREGATION 

A. Discrete description. Master equation 

We consider traffic flow on a single-lane road and study 
the spontaneous formation of a jam regarded as a large 
car vehicle cluster arising on the road. Instead of deal- 
ing with a certain road part of length L and imposing 
some boundary conditions at its entries and exits we ex- 
amine a circular road of length L with N cars moving on 
it. All the cars are assumed to be identical vehicles of 
length Zcar and can make up two phases. One of them is 
the set of "freely" moving cars and the other is a single 
cluster. The cluster is specifled by its size n, the num- 
ber of aggregated cars. Its internal parameters, namely, 
the headway distance /iciust (i-e. the distance between 
the front bumper of a chosen car and the back bumper of 
the following one) and, consequently, the velocity Wciust of 
cars in the cluster are treated as fixed values independent 
of the cluster size n. We note that in the model under 
consideration there can be only one cluster on the road. 
The "free" fiow phase is also specified by the correspond- 
ing headway distance /ifrcc('T-) which, however, already 
depends on the car cluster size n because the larger is 
the cluster, the less is the number {N — n) of the "freely" 
moving cars and, so, the longer is the headway distance 

When a cluster arises on the road its further growth is 
due to the attachment of the "free" cars to its upstream 
boundary, whereas the cars located near its downstream 
boundary accelerating leave it, which decreases the clus- 
ter size. These processes are treated as random changes 
of the cluster size n by ±1 (Fig. ||) and the cluster evolu- 
tion is described in terms of time variations of the prob- 



ability function 'P{n,t) for the cluster to be of size n at 
time t. Then following Mahnke et al. |l^ we write the 
following master equation governing the cluster growth 

dtV{n,t) ^w+{n- l)V{n- l,t) + w^{n + l)r{n + l,t) 
-[w-{n)+w+in)]V{n,t), (1) 

where the cluster size n meets the inequality 1 < n < 
(N — 1) and w+{n) and w^{n) are the transition rates il- 
lustrated in Fig. H and depending generally on the cluster 
size n. The formation and dissolution of the maximum 
possible cluster containing all the cars is described by the 
equation 

dtViN, t) ^w+{N- l)ViN ~l,t)~ {N)V{N, t) , 

(2) 

whereas the emergence of the jam seed, the cluster con- 
sisting of one car called below precluster, obeys the equa- 
tion 



dtV{Q, t) = w^{l)V{l, t) - w+(0)7'(0, t) 



(3) 



Here the function 7^(0, i) is the probability of no clus- 
ter on the road. At the initial time t — Q no cluster is 
assumed to be on the road: 



(4) 



where (5„„' is Kronecker's symbol. The system of equa- 
tions (0)-(^ subject to the initial condition (^) makes 
up the probabilistic description of the cluster formation. 

Special attention should be paid to the question as 
to what the precluster is. The model proposes the fol- 
lowing. When there is no cluster on the road, i.e. all 
the cars move "freely" the velocity of one of them can 
randomly drop down to the velocity Vdnst in the clus- 
ter. Such a car is regarded as the precluster, a size-one 
cluster. When a precluster has arisen its further evolu- 
tion follows the scheme shown in Fig. The preclus- 
ter concept may be justified by recalling the problem we 
deal with initially, i.e. the breakdown processes in mul- 
tilane traffic flow. The cars under consideration actually 
match small vehicle clusters of the synchronized mode, 
macrovehicle, arising in traffic fiow on a multilane high- 
way and comprised of real vehicles moving synchronously 
at different lanes. Therefore the precluster is actually as 
a certain sufficiently small cluster of the synchronized 
mode. Keeping in mind the relatively low threshold of 
the F ^ S transition we will assume the precluster gen- 
eration as well as the precluster dissipation to be inten- 
sive processes so that the "free" flow phase, n = 0, and 
the precluster state, n — \, come into quasi-equilibrium 
on time scales needed for the critical cluster nucleus to 
arise. In particular, in no case the precluster emergence 
limits the cluster evolution, so, the particular details of 
the precluster formation has no substantial effect on the 
traffic breakdown. 

At the next step we should specify the transition rates 
w^in) and w-{n). Let us apply to the optimal velocity 
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FIG. 3: The headway distance ftfroc in the "free" flow phase 
vs the cluster relative volume rj = n/N. A qualitative sketch. 



FIG. 4: The detachment rate 10^(11) vs the cluster size n. A 
qualitative sketch. 



model assuming the velocity v of the "freely" moving cars 
as well as the clustered cars to be determined directly by 
the corresponding headway distance h according to the 
formula 



hP 



(5) 



Here the value -Dopt is the headway distance at which 
drivers feel themselves "free" and their velocity attains 
the maximum Wmax- The parameter p > 1 allows for 
possible forms of the function 'd{h), the greater the value 
of p, the sharper the dependence ■d{h). Case p = 2 is 
often used A car attaches itself to the cluster as 

fast as the distance to the last car in the cluster decreases 
down to the cluster headway /iciust, enabling us to write 
the following ansatz for the attachment rate to the cluster 
of size n > 1 

w+{n) =w^{n) := — . (6) 

«frcc(,"-j — /Iclust 

Applying to a simple geometrical consideration and as- 
suming TV ^ 1 as well as TV — n ^ 1 we get the relation- 
ship (illustrated also in Fig. ^ 



car ~r "clust ) 



(Plim - P) 
p{l - 77) 



(7) 



which together with gives the attachment rate as a 
function of the cluster size n. Here we have introduced 
the following traffic flow parameters: p = N/L being the 
mean value of the car density on the road, its maximum 
possible value pii,„ = l/(?car + ^ciust) for the given road, 
and the relative volume rj — n/N of the cluster with 
respect to the initial volume of the "free" flow state. 

In order to compare the cluster growth due to the car 
attachment with the precluster generation we specify its 
rate in terms of 



w+iO) 



'(0), 



(8) 



where e is a phenomenological factor and we formally set 
n = in expression (^). Keeping in mind the aforesaid 
about the precluster emergence we assume the factor e to 
be about unit, e 1, or at least not to be small enough 



to limit the cluster formation, so its particular numerical 
value is of no importance. 

The rate of the cars escaping from the cluster at its 
downstream front is written as (see also Fig. ||) 



W- (n) 



1 



1 



(n) (j){n) 



T{n) 



To 



(9) 



where the value T(n) can be interpreted as the character- 
istic time needed for the flrst car in the cluster to leave 
it and to go out from its downstream boundary at a dis- 
tance about the headway distance hhcciji) in the current 
"free" flow state. The function (/)(n) allows for the de- 
pendence of the detachment time rin) on the cluster size 
n. We note that expression (|p) is the main original part 
of the model under consideration. 

When the cluster is sufficiently large, n 1, it is rea- 
sonable to regard the characteristic time T(n) « Too as 
a constant (i.e. ^(n) for n ^ 1) as was done in 
papers |l6| for all the values of n. 

For small clusters the r(n)-dependence, however, re- 
quires a special attention. The matter is that the car 
attachment rate w+{n) is considered to be directly de- 
termined by the local characteristics of the "free" flow 
phase and the car cluster. Thus the dependence of the 
attachment rate (n) on the cluster size n arises via the 
headway distance hfj-cc{n) in the "free" flow being a func- 
tion of n, i.e. w+{n) = w+ [/i.free('T-), ^ciust]- Therefore the 
attachment rate is actually an explicit function w+ (p, rj) 
of the mean car density p and the cluster relative vol- 
ume 77 only and, so, exhibits minor variations on scales 
Sn <C N. As will be seen below, exactly this feature is es- 
sential rather than the particular form of w+(p, t]) given 
here for definiteness sake only, because to describe traffic 
breakdown at least one of the kinetic coefficients w+ (n) 
and (n) has to be a direct function of the cluster size 
n for its relatively small values corresponding to the for- 
mation of the cluster critical nucleus. We associate this 
dependence with the escaping rate w-{n) that, in con- 
trast to the attachment rate 11)+ (n), exhibits substantial 
variations in the region n ^ no ~ 20. 

The parameter ng actually divides the car clusters into 
the large cluster group, n ^ rip, for which the escaping 
rate is constant, (j){n) 0, and the group of small clus- 
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ters, n ^ no, whose dissolution is affected substantially 
by the size n. This assumption is based on the fact that 
there should a variety of possible manoeuvres for a driver 
to escape from a sufficiently small cluster on a multilane 
highway when the lanes are not too crowded. 

Expression (H) takes into account this effect via the 
function (/'(n) running from 1 to as the cluster size 
n increases, so, 0(1) ~ 1 and 0(n) ^ as n ^ oo. 
In particular, for a small neighborhood of the precluster 
size, n ~ 1, the value tq of r(n) gives us actually the 
lifetime of the preclusters and is assumed to be less than 
the escaping time from a large cluster, i.e. tq < Too. 
Naturally, for the case of no cluster on the road we have 
to set w-{Q) = 0. The main results will be obtained 
below actually applying to the general properties of the 
dependence w^{n), however for simplicity sake, we will 
adopt the following ansatz for n > 1 



0(n) 



no 



1 



(1 



(10) 



where the exponent g > 1 is regarded as a given constant. 
We point out once more that the dependence of the char- 
acteristic time T(n) on the cluster size is crucial because 
it is responsible for the existence of the metastable "free" 
flow phase. 

The system of equations (|l|)~(P subject to the ini- 
tial condition (^) with the relationships (^), (^), and (||) 
forms the proposed probabilistic model for the car aggre- 
gation. Within this model we will analyze the character- 
istic features of the large cluster emergence and the form 
of the fundamental diagram, i.e. the "flow volume - car 
density" relation in the vicinity of traffic breakdown. In 
particular, in the adopted terms the flow volume j{n) for 
the given traffic flow state, i.e. when a cluster of size n 
arises on the road is written as uB, uM : 



j{n) = (1 - 77)pi? [hhcc[n)\ + r]pi) [/iciust] 



(11) 



Averaging expression ( p, Ij ) with respect to the distribu- 
tion V{n,t) we will get the fundamental diagram j — 
3{P). 



B. Equilibrium distribution 

To clarify the characteristic features of the cluster for- 
mation let us analyze, first, the stationary cluster sizes 
distribution 7^cq(n-)- The system of equations (|l|)-(||) sub- 
ject to the initial condition (^ admits the stationary so- 
lution 7'eq(^) meeting the zero "probability" flux in the 
cluster size space: 

w+ {n - l)Vcqin - 1) - w_ {n)'Pcq{n) = . 

Whence we see that 

Vcqin) _ w+(n~ 1) 



enabling us to write the expression 

Vcqin) cx exp{— f2(n)} , 

where the function U,{n) (called below the car growth 
potential) is specified for n > 2 by the formula 



oo i«-f 



Vt{n) - ^ In [t, 

n' = l 
n I 



w°2(n!)\ 



ro) 



To 



0(n') 



(12) 



Both of the terms in (|lj) vary weakly as the argument 
n changes by one, enabling us to convert sum (^^ into 
an integral with respect to the cluster size n treated as a 
continuous variable: 



fl{n) — Qooin) + f^o('^) 



(13) 



where 



^7oo(^^) 



D,oin) 



dn' In [tooW°^ [/ifrcc(«')]] , (14) 



dn' In 



1 + """^ cj^in') . (15) 

To 



The former term in ( p^ or (|l^), i.e. the component 
f^oo(^) called below the growth potential mainly charac- 
terizes whether a stable car cluster can arise on the road 
under the given conditions and specifies its size because it 
exhibits substantial variations on large scales exceeding 
substantially the size no. By contrast the latter one, the 
component Jig (n) , describes the formation of the critical 
cluster nucleus and, so, the breakdown phenomenon. In- 
deed, as follows from ( [Tc| ) and (^5|) the potential floin) is 
constant for n^ no and, thus, cannot affect the growth 
of a large cluster already formed on the road. Besides, 
within the continuum approximation we have ignored the 
details of the cluster distribution in the region including 
both the points ri = and n = I and expand the cluster 
space n > 1 to the whole axis n > 0. 

Let us, first, analyze the condition of the cluster emer- 
gence. Applying to Fig. ^ we can see that a large cluster 
can arise on the road, in principle, if there exists a value 
of the headway distance he meeting the equality 



TooW°^ [K] = 1 , 



(16) 



Veqin - 1) 



W- (n) 



which will be assumed to hold beforehand. In particular, 
within approximation (^ together with (||) for -Dopt ^ 
/iciust and p = 2 this assumption holds if TooWmax > 2-Dopt 
and the critical headway reads 

'^i^^ = \ (TooWmax + \J (TooWmax)^ " 4L'^pt) , 

whereas for p — > oo and TooWmax > -Dopt we have 

/,(oo) _ 

"'c — ' oo ^max ■ 
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FIG. 5: The attachment rate w'^ [h] vs the headway distance 
h and the stability regions of the "free" flow phase. 



The "free" flow phase will be stable if the initial headway 
distance /ificc(O) > he and unstable otherwise. 

Let us justify these statements. The growth poten- 
tial Q{n) is actually the sum of In [w-{n)/w^{n)] over 
n (see formula (p^)). So, in the region where the inte- 
grand of (|lj) meets the inequality Tocw'^ [h] < 1 and 
the potential Qooin) is an increasing frmction of n, the 
cluster dissolution is more intensive then the car attach- 
ment. Under these conditions the cluster size on the av- 
erage decreases in time. The same concerns the time de- 
pendence of the headway distance hfi-cc{n) in the "free" 
flow phase because the value of hf^cdn) decreases as the 
cluster becomes smaller (Fig. ||), which is also illus- 
trated by arrows in Fig. ||. Since t^oW'^ [h] < 1 for 
h > he any randomly arising cluster tends to disappear 
and, consequently, the "free" flow phase is stable when 
^frcc(O) > he- In this case the potential Uoc{n) possesses 
one minimum located at the boundary point n = (or 
n = 1 what is the same in the continuum description). 

Otherwise, /ifrec(O) < he, there is a region ft.frGc(0) < 
h < he where Tacw'^ [h] > 1 and the car attachment 
rate exceeds that of the cluster dissolution and a cluster 
occurring in the corresponding "free" flow state tends 
to grow, inducing the further increase in the headway 
distance hi^ac{n). In this case the "free" flow phase is 
unstable and the cluster will continue to grow until the 
value of ht-[co{n) reaches the critical point he, where the 
car attachment and the cluster dissolution balance each 
other. Whence it follows, in particular, that the devel- 
oped cluster is of the size ndust obeying the equation 



/jfroc(«clust) — he 



(17) 



and the D.oo{n) has a minimum at the internal point 
n = ficiust- In the present paper we will ignore the exis- 
tence of another region where the equality TqoU'^P' [h] < 1 
also holds for very dense traffic flow, which has been con- 
sidered in papersjl5|, . 

Relationship ^ enables us to rewrite the instability 
conditions in terms of the mean car density p. The criti- 
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FIG. 6: The fundamental diagram in the vicinity of the break- 
down region. 



cal value pc of the car density is the solution of the equa- 
tion /ifrcc(O) — he, whence we immediately get 

^rar "t" ^clust 



Pel = Plii 



^car ~t~ he 



(18) 



Then the stable state of the "free" flow corresponds to the 
inequality p < pei and it loses the stability when p > pei ■ 
In the latter case a large cluster arises on the road whose 
size nciust(p) = f7ciust(p)-^ and relative volume 



Velustip) = 



he + 'car P — Pel 



''clust 



p 



(19) 



In the given analysis we have ignored the dependence 
of the cluster dissolution rate w- [n) on the size n and, 
thereby, the considered picture describes actually the 
"free" flow-cluster transition of the second order. It 
does not allow for the metastable state of the "free" flow 
phase and corresponds to the continuous transition from 
the traffic state of no cluster on the road to the forma- 
tion of a certain cluster whose relative volume changes 
continuously from zero as the car density penetrates 
deeper in the instability region (see formula (|l9[)). Con- 
sequently, this approximation cannot explain the traf- 
fic breakdown and on the phase diagram matches solely 
the stable branches "/" and "c" of the "free" flow and 
the traffic with a developed cluster, respectively (Fig. H). 
Nevertheless, exactly the given approximation describes 
the stable branches of the fundamental diagram and, 
moreover, the metastable branch is a continuation of the 
branch "/" into the instability region. Keeping the latter 
in mind we present also the expression specifying these 
branches: 



j/c(p) 



pd [he + [he + lca.r){Pel ~ P) / P] if P < Pel 
jel -G{p- Pel) I Pel if p > Pel 



where the constants 
3 el = Pci-& \hc\ , 

(he -I- ^car) Pel 



G 



(he — /iclust) Plii 



{pel-d [/ic] - Plim?5 [/iclust]) 
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It should be noted that in obtaining this expression we 
have substituted the maximum probabiUty value ^ciust of 
the cluster size into expression ( p] ) instead of averaging 
it over the distribution Vcq{n). The latter is justified be- 
cause the effect of the cluster size fluctuations is ignorable 
due to iV > 1. 

Now we analyze possible metastable states of the "free" 
flow phase. In order to do this we should take into ac- 
count both the component of the growth potential n{n). 
Since the function ^o{n) exhibits remarkable variations 
in the region n ^ rig only and, thus, the size of the crit- 
ical nucleus also belongs to this region we may confine our 
consideration to clusters whose size n is much less than 
the final cluster size ndust attained after the instability 
development. In addition for the sake of simplicity we 
will regard the value (tqo — tq)/too as a small parameter, 
which enables us to examine solely a small neighborhood 
of the instability boundary, < p — pci <C Pci ■ 

In this case the value of ToqW'^ ['ificc(")] is practically 
constant and can be approximated by the expression 



cluster size, in n/n. 



P- Pel 

I 

Pel 



(20) 



where the coefficient 

(^car H" ^c) 



dlnw°^[/i] 



d\nh 



h—hc 



is about unity, 5 ^ 1, in the general case. In particu- 
lar, for the stepwise dependence 'd[h) (if we set p — 00 
in expression (||)) and Dopt ^ ^car, ^ciust we have the 
rigorous equality g = 1. Expression ( pO[ ) together with 
formula (^^ allows us to represent the dependence of the 
growth potential ri(n) on the cluster size n as 

dVL{n) (Too -To) s {p-pci) 

Pel 

' 5^^^. (21) 



dn To 

_ (Tqc " To) 
To 



no 



no + n 



Pel 



The first term on the right-hand side of (21) is due to 
the increase in the cluster dissolution rate for n ^ no, 
whereas the latter one is proportional to the cluster 
growth rate in the region of large values of n. The re- 
sulting value of the derivative dfl{n)/dn characterizes 
the direction of the cluster evolution. If it is positive, 
dQ(n)/dn > 0, i.e. the potential ri(n) is an increasing 
function of n the cluster dissolution is the dominant pro- 
cess and the cluster of size n tends to dissipate. Other- 
wise, i.e. when dfl{n)/dn < it will grow. 

The former term attains its maximum at n = 0, so, ac- 
cording to ( pT| ) the derivative dfl{n)/dn is negative for all 
the possible values of the cluster size under consideration 
< n ^ ?^ciust when 



P > Pe2 := Pel 



(Too - To) 



1 



5T0 



(22) 



In this case the "free" flow phase becomes absolutely un- 
stable. Under the opposite condition, pd < p < Pe2 there 



I 



a 



o 



o 




cluster size, in n/n„ 



FIG. 7: The form of the cluster growth potential fl(n) in the 
breakdown region, {upper window) the value of —dQ{n)/dn 
proportional to the mean rate of cluster growth vs the cluster 
size n and {lower window) the growth potential r2(n) itself. 
The present figures have been obtained using ansatz (|lO| ) with 
the exponent q = 2 and the chosen value of the vehicle den- 
sity p gives the ratio of the critical nucleus size to the 
characteristic value no equal to Uc/no — 0.8. 



is a certain value Uc at which the derivative d^{n)/dn 
changes the sign (Fig. ^. Setting the left-hand side of 
( ^l| ) equal to zero we get the relationship 



(/)(nc 



Pel 



Pe2 - Pel 

which togather with ansatz ( |l0|) gives the estimate 



(23) 



no 



Pe2 - Pel \ 1 
P- Pel 



- 1 



(24) 



well justified except for small neighborhoods of the 
boundary points pci and pc2- If n < nc the derivative 
is positive and the cluster should decrease in size, i.e. 
the "free" flow phase is stable with respect to the emer- 
gence of such small clusters. However, if a certain cluster 
of size n > has already formed, for example, due to 
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random fluctuations, then it will grow and a large cluster 
of size riciust arises on the road because dQ{n)/dn < in 
the region n > Uc- 

In other words, we have shown that the dependence 
of the dissolution rate w-{n) on the cluster size n ^ hq 
makes the "free" flow phase metastable when the car den- 
sity belongs to the interval p € {pci, Pc2) (branch "m" in 
Fig. ^). The formation of a large cluster, n 3> no, pro- 
ceeds via generation of the critical nucleus whose size Uc 
is estimated by expression (p4|). In order to find the gen- 
eration rate of the critical nuclei and, thus, the break- 
down frequency we should consider the transient pro- 
cesses in the cluster growth, which is the subject of the 
next section. 



C. Continuum approximation. The breakdown 
probability 

In order to apply well developed techniques of the es- 
caping theory (see, e.g., to the analysis of the traf- 
fic breakdown probability we approximate the discrete 
master equation (|l]) by the corresponding Fokker-Planck 
equation. It is feasible because in the case under consid- 
eration the kinetic coefficients w)+(n), u;_(n), first, vary 
smoothly on scales about unity and, second, are approx- 
imately equal to each other, \w^(n) — 'W-{n)\ <C it;_(n). 
The latter conditions enable us to treat the argument 
n as a continuous variable and to expand the functions 
w+{n ± 1), w-{n ± 1), and V{n ± l,t) into the Taylors 
series. In this way and, in addition, taking into account 
expression (^ we reduce equation (|l|) to the following 
Fokker-Planck equation 

T^dtVin, t) = dn [dnV{n, t) + V{n, t)dnn{n)] , (25) 



where the potential Q{n) is given by formula (|13| ) in the 
general form. However, in the case under consideration 
the ratios n/N, [p — Pd)/ Pci, and (tqo — tq)/tq are re- 
garded to be sufficiently small and it is possible to expand 
the potential r2(n) in the tree parameters and to remain 
the leading terms only. In this way we get 




n 
no 



1 - 



2nciust(/o) 



(26) 



where Xc = nc/riQ and we have set J7(0) = 0. Equa- 
tion (|^) transforms into the boundary condition at in- 
finitely distant points that is imposed on the probability 
fiux 

J(n) := -dr,V{n,t)-V{n,t)dnn{n) 

and requires it to be equal to zero, J{oo) — 0. Equa- 
tion (^ describing the precluster generation is reduced, 
in turn, to the zero boundary condition imposed on the 



probability flux J{n) formally at n = 0, i.e. J(0) = 0. 
The latter is justified by the assumed quasi-equilibrium 
between the "free" fiow phase and the preclusters. And, 
finally, the initial condition (^ can be rewritten as 



dnV{n,t) = 1. 



When the car density belongs to the interval p S 
{Pci,Pc2) and the "free" flow phase is metastable the 
form of the growth potential 0(n) in the region n uq is 
shown in Fig. ^ The the "free" flow phase being in quasi- 
equilibrium with preclusters matches the local minimum 
at n = separated from the region of the stable cluster 
growth n > Uc hy the potential barrier flc- The value of 
this the potential barrier is estimated as 



where the function 



To 



-71o ^ 



(27) 



dx ; 



dx 



(28) 



In particular, for ansatz (^^ with the exponent g = 2 it 
expression (U% becomes 



fir 



(Too - To) 



Xr 



-no- 



To ' (1 + Xc)^ 

moreover, in the limit <C 1 we have 
uj[xc] — 2^^c I where the constant r — - 



(29) 



d(l)[x] 



dx 



as follows from expression (p8D, and the general formula 
for the potential flf. can be written as 



(To, 



2ro 



To) 2 

rriQX^ . 



(30) 



In the same limit expression (^) gives us 
(Pc2 - P) 



r{pc2 - Pel) 



(31) 



The main much deeper minimum of the potential fl(n) 
is located a.t n — rtciust ^ "o- 

We have demonstrated that a precluster must climb 
over the potential barrier Qc at the point ric to convert 
into a large stable cluster. It is implemented through 
random fluctuations carrying the cluster size up to the 
critical value ric- In this terms the traffic breakdown 
is the classical escaping from a potential well described 
by the Fokker-Planck equation (p5|). The latter analogy 
enables us to write down the estimate for the frequency 
Vhd of the traffic breakdown processes depending on the 
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given vehicle density in the "free" flow state. Namely, as 
shown in Appendix 



1 



3 

Tq\2 



V27rnoroo ^ To 



X exp 



(Too - To) 



To 



nouj[xc] } , (32) 



which is well justified for the car density p belonging 
to the interval pd < p < Pc2 except for a certain suf- 
ficiently small neighborhoods on the critical points pd, 
Pc2- Ansatz (10) with the exponent q = 2 together with 
formula (23) enables us to rewrite expression (p^) as 



1 




Too-T0A2^^_^^^| 



(Too - To) (1 - A)2 
Uq- 



To (1 + A^) 

Here we have introduced the quantity 



A 



P- Pel 
Pc2 - Pel 



(33) 



(34) 



treated as a dimensionless overcriticality measure show- 
ing how deep the system penetrates into the metastability 
region, A = corresponds to the value pci of the vehicle 
density where a jam can emerge in principle and A = 1 
matches the vehicle density pc2 after exceeding which no 
traffic states except for jams can exist at all (Fig. 



D. Frequency of traffic breakdown during a fixed 
time interval 



Experimentally traffic breakdown is typically analyzed 
detecting a significant drop in the vehicle speed during a 
certain fixed time interval Tobs about several minutes and 
then drawing the relative frequency of these events vs the 
traffic volume |, |l^, |l^, |l|. In order to compare this 
representation with the obtained results let us consider 
them in more details. 

As follows from expression ( ^ ) the density interval 
(pci, Pc2) inside which the traffic jam emerges by the nu- 
cleation mechanism is of the thickness 



{Pc2 - Pel) = p 



cl- 



(T qo - To) 
gTQ 



According to the experimental data ^ |T^, |l^, |l3| the 
thickness of the traffic volume interval inside which the 
traffic breakdown demonstrates the probabilistic behav- 
ior is about its low boundary in magnitude. So we have 
to regard the ratio (tqo — to)/to also as a value about 
unity: 



(tqo - To) 

5To 



1. 



(35) 



Thereby, setting uq = 20 we get the conclusion that 
in the general case where ~ uq the potential bar- 
rier ~ 5 corresponding to the exponential factor 
exp{— ilc} ~ 0.7 X 10~^. Then setting Tqo ~ 2 sec and 
estimating the preceding cofactor as 1/ (•\/27rnoToo) we 
find the characteristic rate of the traffic breakdown be- 
ing about 1/50 min~^ in the general case. So the real 
traffic breakdown events seem to be observed in cases 
where the vehicle density comes to the upper boundary 
Pc2- The latter allows us to confine our analysis formally 
to the limit case 

< 1 {pc2 - p) < {Pc2 - Pel)- (36) 

Then estimating the probability J^bd of detecting a 
traffic breakdown during the observation time interval 
Tobs as JFbd = Tobs'^bd wc obtain from (^2|) the expres- 
sion 



^bd(A) 



Tbd 



■2ni exp{-f7c} 



where 



2rTo 



and we have introduced the time scale 



Tbd 



2^To 



f(To 



:»^oTo< 



(37) 



(38) 



(39) 



giving us the characteristic time of the breakdown emer- 
gence. In deriving ( p7| ) we have also taken into account 
formulae ([lo| ) , (|30| ) , and (|3^) and remained directly num- 
ber 2 as cofactor because the maximum of the function 
z^/^ exp(— z) is about 0.43. Naturally, we have to confine 
ourselves to such values of the vehicle density for which 
^hd < 1 because traffic flow with higher values of the 
vehicle density cannot exist on these time scales. For 
the following values r = 2, the ratio (too — to)/to ~ 1, 
Too ~ 2 sec, and uq ^ 20 expression ( p9| ) gives us the es- 
timate Tbd 1 min of the characteristic breakdown time. 
It should be pointed out that the latter estimate does not 
contradict the evaluation of the breakdown rate given at 
the beginning of the present section because it holds only 
in the region ilc ^ 1. 

Figure || illustrates the obtained results depicting the 
breakdown probability for different values of the obser- 
vation time Tobs measured in units of Tbd vs the depth 
of penetrating into the metastability region. It should 
be pointed out that in drawing Fig. g| we have applied to 
formula ( ^ ) rather than (^) in order to have a possibil- 
ity to go out of the frameworks of the formal limit (|3^). 
The latter is considered hear to clarify the obtained re- 
sults only. To make the form of the jrbd(A)-dependence 
more evident we apply again to the formal limit case ( |3^ ) 
assuming the ratio m := Tobs/Tbd to be a large param- 
eter. Than analyzing a small neighborhood of the point 
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ansatz ( p^ ) with the exponent q = 2 




overcriticaUty measure, (p - P^^lip^^ ~ 1^^) 



volume from about 1500 veh/h/1 (vehicles per hour per 
lane) up to 3000 veh/h/1. The real dynamics of traffic 
breakdown near bottlenecks and the developed structure 
of the congested traffic flow are sufficiently complex as 
it was exhibited by Kerner In particular, Kerner 

demonstrated that the synchronized mode of traffic flow 
in the vicinity of highway bottleneck is locally metastable 
under the discharged downstream traffic flow of volume 
j varying in the same interval. The latter enables us to 
estimate the detachment time Tqo playing the significant 
role in the presented model. In fact ignoring the velocity 
'^(/iciust) of cars in the cluster as well as the headway dis- 
tance /iciust in ansatz ^ we get that the lower boundary 
Pel of the metastability region meets the traffic volume 



[he] = — - 1800 veh/h/1. 



FIG. 8: The traffic breakdown probability vs the depth (p — 
pci)/{pc2 — Pel) of penetration into the metastability region. 



p^ specified by the equality 



A* 



2r ln(2m)To 



(too - To)no 



we get 



^hdip) = exp 



P - Pn 



(40) 



(41) 



where the vehicle density scale 



Pel) 



rro 



21n(2TO)(roo - 7-0)710 



(42) 



Therefore, in a rough approximation the J-hd{p)- 
dependence is a simple exponential function whose scale 
p„ is approximately a constant value (because the func- 
tion ln(m) shows week variations for m ^ 1). Changing 
the observation duration Tots practically shift the cut-off 
point p^ only. 



III. COMPARISON WITH EXPERIMENTAL 
RESULTS 

Experimental investigations of the traffic breakdown 
regarded as a probabilistic phenomenon have been car- 
ried out by several authors (see, e.g. Ref. ^, |l^, |l^, 
|l2|, |l^). Elefteriadou et al. Q actually pointed out 
the fact that traffic breakdown at ramp merge junctions 
occurs randomly without precise relation to a certain 
fixed value of traffic volume. A more detailed analysis 
of the breakdown probability has been fulfilled in pa- 
pers |l^, ^ |l^. These observations show that traf- 
fic breakdown can occur inside a wide interval of traffic 



Whence we find immediately the estimate Tqo ~ 2 sec, 
which is i n ag reement with the value adopted previously 
in papers |l^ |l^. In this section, as it has been done in 
the previous one, we use the estimates of the quantities 
no ~ 20 according to the experimental data depicted in 
Fig.|l|, set (too — To)/ro ^ 1 from the general considera- 
tion. 

In order to compare the obtained results and the avail- 
able experimental data we have applied to the latest ma- 
terials presented in detail by Lorenz & Elefteriadou [p^ . 
The breakdown phenomenon was investigated in traffic 
flow near two bottlenecks of Highway 401, one of the pri- 
mary Toronto traffic arteries. The detectors were located 
right after the on-ramps within several hundred meters 
downstream. So the dynamics of traffic breakdown ob- 
served at these places seems to be mainly due to local 
internal properties of traffic flow discussed in the present 
paper. The complex spatial structure of the induced con- 
gested phase including moving wide and narrow jams re- 
ported by Kerner |]l3| should emerge above the detectors 
upstream. We consider in detail the data obtained for 
one these bottleneck (site "A" in paper |l^ ) . The paired 
detectors were located in each of the tree lanes and were 
instrumented to provide vehicle count and speed esti- 
mates continuously at 20-second intervals. A breakdown 
event was fixed via the velocity drop below 90 km/h, the 
middle point of a certain gap in the velocity field visu- 
ally separating the congested and free traffic flow states. 
Besides, only those disturbances that caused the aver- 
age speed over all the lanes to drop below 90 km/h for 
a period of five minutes or more were considered a true 
breakdown. The latter enabled the authors to filter out 
large amplitude fluctuations in the mean vehicle veloc- 
ity not leading to the traffic breakdown. Fig. jo] exhibits 
the obtained probability (relative frequency) of the traffic 
breakdown events during 5-minute and 15-minute inter- 
vals vs the traffic volume partitioned within 100 veh/h/1 
steps. We note that Fig. g does not show the available 
1 minute interval data because the corresponding break- 
down probability is not significant for all the observed 
values of traffic volume except for the upper boundary 
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FIG. 9: The traffic breakdown probability during 5-minute 
and 15-minute intervals vs the traffic volume (after Lorenz & 
Elefteriadou [Q). The continuous curves present the depen- 
dence (|3^) fitting these experimental data in the frameworks 
of the replacement A <^ (j — jci)/{jc2 — jci) (see expres- 
sion (^) ) under the following values of the critical traffic vol- 
umes jci = 1200 veh/h/1 and jc2 = 3400 veh/h/1, the char- 
acteristic breakdown time Tbd ~ 2.5 min, and the product 
no(roo — to)/to = 25. 



2800 veh/h/1, that can be due to its rare occurrence. 

The continuous curves in Fig. |9| present our attempts 
to fit the obtained theoretical dependence in the given 
experimental data. Namely these curves describe the 
breakdown probability estimated as J^bd = Tobs^'bd, 
where the latter cofactor is given by the expression (^) 
within the replacement A <;= (j — jci)/{jc2 — jci) (see 
expression (p3)) and we have used the following values 
of the critical traffic volumes jci = 1200 veh/h/1 and 
jc2 = 3400 veh/h/1, the characteristic breakdown time 
Tbd = 2.5 min (see expression ^), and set the product 
''t-o{too — T())/tq = 25. Keeping in mind the aforesaid 
all the adopted values are quite reasonable including the 
estimate jc2 — 3400 veh/h/1 seeming at first glance ex- 
tremely high. Indeed, this value is no more than a result 
of approximating the j(yo)-dependence by a linear func- 
tion formally into the region of high vehicle densities. 



IV. SUMMARY AND CONCLUSION 

We have considered the traffic breakdown phenomenon 
regarded as a random process developing via the nucle- 
ation mechanism. The origin of critical jam nuclei pro- 
ceeds in a metastable phase of traffic flow and seems to 
be located inside a not too large region on a highway, 
for example, in the close proximity of a highway bottle- 
neck H, ^. The induced complex structure of the con- 
gested traffic phase is located upstream the bottleneck 
p3| . Keeping these properties in mind we have applied 
to the probabilistic model regarding the jam emergence 
as the development of a large car cluster on highway. In 



these terms the traffic breakdown proceeds through the 
formation of a certain car of critical size in the metastable 
vehicle flow, which enabled us to confine ourselves to the 
single cluster model. 

We assumed that, first, the growth of the car cluster 
is governed by attachment of cars to the cluster whose 
rate is mainly determined by the mean headway distance 
between the cars in the vehicle fiow and, may be, also by 
the headway distance in the cluster. Second, the clus- 
ter dissolution is determined by the car escape from the 
cluster whose rate depends on the cluster size directly. 
To justify the latter assumption we apply to the mod- 
ern notion of the traffic flow structure (see Ref . |l|, |^, |^) . 
Namely, the jam emergence goes mainly through the se- 
quence of two phase transitions: free flow —^ synchronized 
mode — > stop-and-go pattern 0. Both of these transi- 
tions are of the first order, i.e. they exhibit breakdown, 
hysteresis, and nucleation effects Therefore consid- 
ering the final stage of the jam emergence we have to 
regard the synchronized mode as the metastable phase 
exactly inside which a critical jam nucleus appears due 
to random fluctuations. The synchronized mode is char- 
acterized by strong multilane correlations in the car mo- 
tion and, as a result, all the vehicles in a certain affective 
cluster spanning over all the highway lanes move as a 
whole. So the proposed probabilistic description deals 
with actually macrovehicles comprising many individual 
cars. The available single-vehicle experimental data |^ 
present the correlation characteristics of the synchronized 
mode which have enabled us to estimate the character- 
istic dimension uq ^ 20-30 of the car cluster entering 
the dependence of the car detachment rate on the cluster 
size. Namely for small car clusters, n ^ uq, the charac- 
teristic detachment time tq should be substantially less 
than this time Too for large clusters, n ^ uq. 

We have written the appropriate master equation for 
the cluster distribution function and analyze the forma- 
tion of the critical car cluster due to the climb over a cer- 
tain potential barrier. The inequality no 3> 1 has opened 
us the way to convert from the discrete master equation 
to the appropriate Fokker-Plank equation and find all the 
required characteristics of the traffic breakdown. 

The obtained results are compared with the available 
experimental data and, in detail, with the probability of 
traffic breakdown in the vicinity of bottlenecks vs the 
traffic volume presented by Lorenz & Elefteriadou p2[ . 
It turned out that the theoretical curves can be fitted 
closely to the given experimental data using values of the 
main parameters chosen based on the general properties 
of the traffic fiow not related directly to the breakdown 
dynamics. In particular, first, we have demonstrated that 
the characteristic internal time scale Tbd of the break- 
down development is about Tbd ~ it-qToo (we recall that 
Too ~ 2 sec is the characteristic time during which a car 
can individually leave a cluster) . Whence we get the esti- 
mate of the breakdown timescale about one minute. The 
latter justifies the widely used probabilistic technique of 
the breakdown investigation based on fixing this event 
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during a time interval of several minutes. Second, the 
proposed model explains why the trafhc breakdown as a 
probabilistic phenomenon is observed inside a sufficiently 
wide interval of the traffic volume, namely, the thickness 
Aj of this layer can attain its low boundary jd in mag- 
nitude. The matter is that Aj/jci ^ (too — to)/to ^ 1. 

Concluding the aforesaid we state that traffic break- 
down is a mesoscopic process, as it must be for the syn- 
chronized mode, whose characteristic spatial and tem- 
poral scales correspond to car clusters made of a large 
number of vehicles. 
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APPENDIX A: ESCAPING RATE FROM A 
BOUNDARY WELL 



In section |II C| we have obtained the Fokker-Plank 
equation ( p5| ) governing the evolution of car clusters 
treated as random wandering in the space of their size 
n. It has turned out that near the threshold the preclus- 
ter domain is separated from the large cluster region by 
a potential barrier 17^, so, the formation of the clustered 
phase should proceeds through the nucleation mechanism 
(Fig. 1^). In other words, for a large cluster to emerge 
on the road its critical nucleus ric has to arise via ran- 
dom fluctuations of the cluster size in the precluster re- 
gion. Thereby in order to describe the cluster formation 
we need the expression specifying the rate of the criti- 
cal nucleus generation, being the subject of the present 
appendix. 

Mathematically the description of the critical nucleus 
generation is equivalent to the problem of a particle es- 
caping form the corresponding potential barrier (Fig. 10). 
Thereby the rate of the critical nucleus generation, i.e. 
the frequency of the traffic breakdown z^bd is represented 
in terms of the probability density T{t") for this particle 
to escape from the potential well at a given time t pro- 
vided initially, t = 0, it has been placed near the local 
minimum (here n = 0). Namely 



v^,A = ^(+0) , 



(Al) 



where the value -1-0 of the argument t means that we con- 
sider time scales exceeding substantially the duration of 
all the transient processes during which the distribution 
of the particle inside the potential well attains locally 
quasi-equilibrium. 

Since the potential relief under consideration is rather 
special we prefer to recall briefly the way of deriving the 
probability Tif) referring a reader to the specific litera- 
ture (see, e.g., book |l^) for details. 

The concept of potential well implies that the barrier 
is sufficiently high, fie ^ 1, therefore the particle can 



climb over it due to rare fluctuations lifting the particle 
to points at the potential barrier where 0(?i) ^1. If 
such an event does not lead to escape the particle will 
drift back to the the neighborhood of the local mini- 
mum n = whose thickness is specified by the inequality 
rj(n) < 1. Thereby the subsequent attempts of escap- 
ing may be considered as being mutually independent. 
After the particle has climbed over the barrier the force 
—d£l(n) I dn carries it away to distant points, making the 
return impossible. So from this point of view we may 
refer to the particle being inside the potential well or 
having escaped from it as two its possible states without 
specifying the particular position. Therefore the proba- 
bility Vit — t') that the particle remains inside the well 
at time t if it has being placed in it at time t' obeys the 
equation: 

■p{t) =V{t~t')V{t') forO<t'<t. 
Whence we get the general expression for the function 

m 

t 



V{t) = exp 

^ Tlifc 

where Tiife is a certain constant specified by the particular 
properties of a potential well. The latter formula gives us 
immediately the general form of the escape probability 



dVit) 1 / t 

m = ^ = cxp 

at -Hifc V Tiifc 



(A2) 



In order to find the lifetime we will deal with the 
Laplace transform J-l{s) of the escape probability J-{t) 



oo 

T^U) := / dt exp(-st)T{t) = 

J 1 + s 



STlifo 



(A3) 



whence it follows that in the expansion of J-l{s) with 
respect to s around the point s = 



J^l(s) 1 - STlif, 



(A4) 



the first order term directly contains the desired lifetime 
as the coefficient. 
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Following the standard approach |Q we reduce the 
escaping problem to finding the first passage time prob- 
ability. In other words, we assume the particle never to 
come back to the potential well if it after climbing the 
barrier reaches points where flc — n{n) > 1 (Fig. 10). 
The particle may be withdrawn from the consideration 
or, what is the same, it will be trapped when reaches for 
the first time any fixed point n* in this region. The time 
it takes for the particle to reach the point n* after over- 
coming the barrier at the critical point is ignorable in 
comparison with the characteristic wait of critical fluctu- 
ations. Thereby the function !F{t) specifies actually the 
probability of passing (reaching) the point n* for the first 
time at the time moment t. This construction enables us 
to introduce a more detailed relative function J^{n, t) giv- 
ing the probability for the particle initially placed at the 
point < n < n* to reach first the right boundary n* 
of the region under consideration at the time moment t. 
The left boundary n = is impermeable for the particle. 
Then using the standard technique based on the back- 
ward Fokker-Planck equation conjugated with Eq. (l2~ 
we the governing equation for the function J^l(ji, s 



t^sTl = OITl - [dnn{n)] [dn^L] (A5) 
subject to the boundary conditions 

TL{0,s)^TL{n*,s) = l. (A6) 

Whence it directly follows that the first order term Lp{n) 
in the expansion of the Laplace transform J-'L{n, s) with 
respect to s 



T{n, s) = 1 — SLp{n) , 
obeys in turn the equation 



subject to the boundary conditions 

a„^(0) = and Lp{n*)=Q. (A8) 
The solution of the system ( |A7| ) and ( A8) has the form 



^(n) = Too / dn' e"("') j dn" e'^f"") (A9) 



and the value (p{0) gives us the desired lifetime: 



nifc = ^{0) . 



(AlO) 



Inside the potential well the function if{n) takes practi- 
cally a constant value mainly contributed by the points 
n" belonging to the well bottom, i.e. to the region 
ri(n") ^ 1 and by the points n' located near the top 
of the potential barrier where ^l{nc) — ri(n') < 1. This 
feature leads us immediately to the approximation 



(All) 

which is the main result of the present appendix. 

In particular, for the potent ial Q (n) specified by ex- 
pression (^ ) or (|2^)) formula ( All ) gives 



To 



Too - To 



X (l-0M)-i(|0'[a;,]|)~2ef^(»=). (A12) 



(A7) 



Formulae (|Al|), (|A2[ ), (|A12|) , and (|27|) give us expres- 
sion (11). 
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